Growth-rate dependence of a bacterial genetic oscillator dynamics 
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Gene networks exhibiting oscillatory dynamics are widespread in biology. The minimal regulatory 
designs giving rise to oscillations have been implemented synthetically and studied by mathemati- 
cal modeling. However, most of the available analyses generally neglect the coupling of regulatory 
circuits with the cellular "chassis" in which the circuits are embedded. For example, the cell macro- 
molecular composition of fast-growing bacteria changes with growth rate. As a consequence, impor- 
tant parameters of gene expression, such as ribosome concentration or cell volume, are growth-rate 
dependent, ultimately coupling the dynamics of genetic circuits with cell physiology. This work 
addresses the effects of growth rate on the dynamics of a paradigmatic example of genetic oscillator, 
the repressilator. Making use of empirical growth-rate dependences of parameters in bacteria, we 
show that the repressilator dynamics can switch between oscillations and convergence to a fixed 
point depending on the cellular state of growth, and thus on the amount and quality of the food it 
is fed. The substrate of the circuit (type of plasmid or gene positions on the chromosome) also plays 
an important role in determining the oscillation stability and the growth-rate dependence of period 
and amplitude. This analysis has potential application in the field of synthetic biology, and suggests 
that the coupling between endogenous genetic oscillators and cell physiology can have substantial 
consequences on their functionality. 
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I. INTRODUCTION 

Oscillatory behavior is widespread and fundamentally 
important in biological systems, from circadian clocks to 
cell-cycle control [l[ . At the level of a single cell, oscilla- 
tions can be sustained by regulatory circuits whose basic 
components are nucleic acids, proteins, and biochemical 
interactions 0, However, naturally occurring regu- 
latory systems are complex. Thus, our path to ratio- 
nalize their behavior has to pass through simplifications 
and basic underlying principles The study of mini- 
mal circuits, by mathematical modeling and experimen- 
tal synthetic engineering, allows a quantitative approach 
to gene expression from regulatory circuits • 

This simplifying approach is not only useful for study- 
ing the extant regulatory circuits. Understanding the ba- 
sic principles and dynamical properties of regulatory net- 
works makes it possible to design and produce synthetic 
circuits that can perform specific functions in a pre- 
dictable manner j9l4l2l|. with potentially relevant future 
biotechnological and biomedical applications [ll|i [13 ■ ^ 
paradigmatic example of a synthetic genetic oscillator, 
and one of the earliest in vivo realizations [3] , is the so- 
called "repressilator", a three-gene cyclic circuit where 
each gene protein-product represses the synthesis of its 
successor (Fig. [T]A). 

From a physics standpoint, dealing with few well- 
characterized isolated systems is necessary for building 
testable models. Unfortunately, every genetic circuit, en- 
dogenous or synthetically implemented in a living cell. 
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cannot be truly considered isolated from cellular pro- 
cesses. These processes are strongly affected by the phys- 
iological state of the cell [Hj, and thus by the environ- 
mental conditions and type/availability of nutrients. For 
example, the macromolecular composition of bacterial 
cells in steady-state exponential growth changes substan- 
tially with growth rate (the rate of cell proliferation) [l6l - 
Importantly, in this case, growth rate appears 
to recapitulate most of the physiological changes. The 
growth-rate dependence of important global cell param- 
eters related to gene expression (such as ribosome and 
RNA polymerases concentrations), couples the dynam- 
ics of gene expression with the cell state, with relevant 
consequences for genetic circuit functioning [l^ [l^ 2^ ■ 

For E. coll, the best studied bacterial species, the 
growth rate dependences of various cellular parameters 
were evaluated quantitatively in a study by Klumpp and 
coworkers [T5| . compiling results from multiple experi- 
ments. Leveraging on these empirical data, they ana- 
lyzed the growth-rate dependence of the steady state of a 
constitutively expressed (i.e. unregulated) gene and few 
other simple genetic circuits, such as the self-regulator 
and the toggle switch. Their modeling strategy will be 
briefly reviewed in section III Al Beyond steady state, 
knowing the empirical growth-rate dependence of gene 
expression parameters allows in principle to revisit the 
dynamical properties of genetic circuits, introducing the 
physiological cell state as a new player. 

This work addresses for the first time the growth-rate 
dependence of a biological oscillator dynamics, and con- 
siders the effect of cell state on a repressilator, integrated 
on either a plasmid or the E. coli chromosome. After 
reviewing the conditions for oscillations and the repressi- 
lator dynamics at fixed-growth rate (and extending some 
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known results), wc will show how the physiology of the 
cell can alter qualitatively and quantitatively the dynam- 
ical properties of the system. The main result is that the 
conditions for stable oscillations, as well as their ampli- 
tude and period are growth-rate dependent. This implies 
that a genetic oscillator can display distinct dynamical 
behaviors in different environments and nutrient condi- 
tions. Additionally, the circuit substrate, e.g. type of 
plasmid or chromosomal position of the genes, also con- 
tributes to its dynamics, through gene dosage. Specifi- 
cally, the range of growth rates in which oscillations are 
observable will vary with substrate, in a manner that 
could be exploited both biologically and technologically. 
Our predictions, although based on a simplified model, 
are experimentally testable in a straightforward way. 

The paper is organized as follows. Section |TT] reviews 
the modeling strategy adopted to include the growth rate 
as a variable in gene expression (Subsection III A) ) and ex- 
plains how this strategy can be applied to a mathemati- 
cal description of the repressilator (Subsection lllBI). Ap- 
pendix Rl reviews the Cooper-Helmstetter model [21[, an 
empirical model of DNA replication in fast-growing bac- 
teria that is an essential integral part of our approach, 
and justifies in more detail some of the model assump- 
tions presented in Section |lll Section IIIII contains the 
quantitative analysis of the circuit, and in particular of 
the role played by the cellular state of growth on the 
dynamics. More specifically. Subsection IIII Al analyzes 
the symmetric repressilator, showing the possibility of a 
growth-mediated switch between oscillations and conver- 
gence to a stable fixed point. Subsection IIIIBI focuses 
on the asymmetric case, with particular emphasis on the 
role of gene chromosomal position on the circuit behav- 
ior. Finally, the last section discusses the implications of 
the results from both a biological and a physical modeling 
standpoint. 



II. MODEL 

A. Growth-rate dependence of gene expression 

This section reviews the approach of Klumpp and 
coworkers [isj . and the necessary assumptions in order 
to extend it to the analysis of circuit dynamics. This 
modeling strategy constitutes the basis for our descrip- 
tion of the repressilator, introduced in section HTBI 

In absence of regulation, the dynamics of messenger 
RNA (niRNA) levels and protein concentrations, denoted 
with i and p respectively, is described by two equations 
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where ai and ap are the transcription and translation 
rates respectively, and I3p the degradation rates of niR- 
NAs and proteins, and V the cell volume. The parameter 



g is the gene copy number. If the gene substrate is a plas- 
mid, g represents the (mean) plasmid copy number. For 
genes integrated on the chromosome, the copy number 
can vary because of DNA replication, especially at fast 
growth rates where multiple copies of the genome are 
replicated at the same time. This phenomenon whereby 
gene dosage is modified by DNA replication is described 
by the Cooper-Helmstetter model [2l|, and reviewed in 
Appendix lA 11 Note that for fast E. coli growth, g in- 
creases with decreasing distance m from the replication 
origin (illustrated in Fig.[Tj3), since the cell engages over- 
lapping rounds of DNA replication in order to allow fast 
growth |22| . 

In bacteria, proteins are typically stable, with a life- 
time longer than the cell cycle, while mRNAs have a 
lifetime of just a few minutes [23. Therefore, the loss 
of protein is mainly due to dilution through growth 
and cell division, so that an effective degradation rate 
/3p = /i ln2 (where is the growth rate) can be safely 
used in most cases. On the other hand, the fast time- 
scale of mRNA dynamics allows a quasi-equilibrium ap- 
proximation. Thus, Eq. [T]can be reduced to 
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Rescaling time with the dilution rate /3p, all the growth- 
rate dependence can be factorizcd in a single term F{fi), 
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Normalizing the growth-rate dependence F(n) such 
that i^(/i) = 1 for /i = 1 db/hr (i.e. doublings per hour), 
the parameter p* represents the steady-state concentra- 
tion of the unregulated protein at 11= 1 db/hr. Thus, it 
is a measure of the degree of basal expression. Follow- 
ing [isj , we will refer to it as the "promoter strength" , 
although this terminology might be slightly misleading, 
as the term actually includes non-transcriptional param- 
eters, such as the translation efficiency or the gene copy 
number at /i = 1 db/hr. 

The function -F(/i) is fully determined by the empirical 
dependences on growth rate of different cellular param- 
eters. For E. coli, the numerical values of these pa- 
rameters were collected for five growth rates ^ between 
0.6 db/hr and 2.5 db/hr (see Table SI of [H). Following 
these values, F{^) decreases in a weakly nonlinear fash- 
ion in the available range (interpolation will be used in 
the following when needed). 

The empirical growth-rate dependences that define 
F{^) arc based on experimental measurements of the 
average cellular properties in a growing cell population. 
Thus, they are in principle also dependent on the age dis- 
tribution (where "age" stands for stage of the cell cycle) 
of the population. However, the differences between aver- 
ages calculated over a cell cycle and over a cell population 
with the age distribution determined by the exponential 
growth are not quantitatively significant, as discussed in 
detail in Appendix |^ Moreover, wc will not consider 
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FIG. 1: (Color online) A. Scheme of the repressilator. 

Each node represents the expression level of one gene, links 
represent repressive regulatory interactions, as described by 
Eq. B. Scheme of the E. coli chromosome, illus- 

trating the coordinate m, used for the positions where a gene 
can be inserted. 



explicitly the cell-cycle dynamics of cellular parameters, 
which, for sufficiently long time scales, can be averaged 
out using for example effective gene dosage or effective 
protein dilution rate. As a consequence, Eq. |3] can ac- 
curately describe the dynamics of protein concentration 
on time scales longer than the generation time [1^ , as it 
will always be the case in the following. 



B. The repressilator 

The repressilator circuit (Fig. [T]A_) can be described by 
three equations, one for each gene, based on the gene 
expression model discussed above: 

xi = x\Fi{^)R{x2/ki) ~ xi 

X2 = xlF2{fi)R{x3/k2) ~ X2 

X3 = xlF3{fl)R{xi/k3) -X3 . (4) 

Transcriptional repression is modeled phenomenologi- 
cally using a nonlinear function of the repressor concen- 
tration (Hill function) [5| |6!| ; 



where the Hill coefficient n defines the degree of coopcr- 
ativity (determining the steepness of R), while the dis- 
sociation constant k specifies the repressor concentration 
at which the production rate is half of its constitutive 
value. 

Assuming similar dissociation constants for the three 
promoters fci ~ ^2 — fcs = fc, the protein concentrations 
can be rescaled with k 



xi = xlFi{fi)R{x2) ~ xi 

X2 = xlF2{fJ.)R{x3) - X2 

X3 = xlF3{n)R{xi) - X3. 



(6) 



With these notations, the promoter strengths x* have 
dimensionless units, as they represent the steady-state 
concentration of protein i in absence of repression at /i = 
1 db/hr rescaled with the dissociation constant k. The 
growth rate functions Fj(/z) can be gene specific due to 
the dependence of gene dosage gt on the chromosomal 
gene position. 



III. RESULTS 

A. Growth rate affects qualitatively and 
quantitatively the dynamics of a symmetric 
repressilator 

We start the analysis from the simplified case of a 
symmetric repressilator, i.e. the three genes have ap- 
proximately the same production/degradation rates for 
mRNA and proteins, as well as the same growth-rate 
dependence of parameters (this is attained for example 
when all genes are integrated on a single plasmid or at a 
similar distance from the origin of replication on a chro- 
mosome). This is the most studied case, and was ap - 
proached with diflEerent mathematical descriptions [1, [25l - 
|27| . The symmetric approximation was originally pro- 
posed to explain the behavior of the synthetic realization 
in vivo 

In the symmetric case, the functions encoding 
the growth-rate dependence F{fi) and the promoter 
strengths x* are exactly the same for each of the three 
genes, simplifying Eq. [6]to 



Xi = x*F{^i)R{xi+i) - Xi , 
where i £ [1,3] and Xi = xi. 



(7) 
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1. At fixed growth rate, the dynamics is determined by 
promoter strength and cooperativity 

Before addressing the efi^ects of the growth-rate depen- 
dence of parameters, we characterize the circuit dynamics 
at fixed growth conditions, for simplicity at fj, ^ 1 db/hr 
where F{ii) = 1. In general, the symmetric repressilator 
can display stable oscillations, arising through a Hopf 
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cooperativity n time [units of |3p] 



FIG. 2: (Color online) Stability diagram of the symmetric repressilator. The left panel shows the parameter space 
giving rise to oscillations. The promoter strength x* is dimensionless, measured as protein concentration at /i = 1 db/hr in 
absence of repression in units of the dissociation constant k. High repression cooperativity and strong promoters can ensure 
stable oscillation. The right panel plots are illustrative examples of the dynamics in correspondence with two sets of parameter 
values ({n = 4,2:* = 1.8} and {n = 4,2;* = 5}), showing the time evolution of the three protein concentrations 21,22,2:3 (in 
units of k). The red continuous curve describes the dynamics of 21, the blue dashed curve 22, and the orange dot-dashed curve 
23 (same colors of the scheme in Fig. [HA-). 



bifurcation [2^. More specifically, Eqs. [7]have an oscil- 
latory solution if the condition for cooperativity n > 2 
is satisfied, as can be shown in a strai ght forward way 
considering the symmetry of the system [2g . 

For a given steepness of the repression function satis- 
fying the condition n > 2, the stability of the limit cycle 
is solely determined by the promoter strength x* . The 
values of x* that ensure stability of the oscillatory state 
can be calculated using linear stability analysis and the 
resulting stability diagram is shown in Fig. [2] Essentially, 
an increase of either cooperativity or promoter strength 
can help the stabilization of oscillations. This result is 
common to all the different repressilator descriptions pro- 
posed in the literature 

Furthermore, the two parameters n and x* determine 
the amplitude and period of oscillation (Fig. |3|). The ap- 
proximately linear and logarithmic dependences of ampli- 
tude and period respectively that emerge from numerical 
integration of Eqs. [7] can be rationalized by the follow- 
ing rough but conceptually simple argument. Each gene 
in the repressilator tends to oscillate between two states 
corresponding to its maximally repressed state and its 
maximally activated state. In a simplified situation of a 
step repression function (n ^ 1) that switches the tar- 
get gene on and off after its equilibration, the protein 
concentration would go from Xi ~ (fully repressed) to 
Xi ~ X* (fully activated), thus leading to a linear depen- 
dence on promoter strength of the oscillation amplitude. 
On the other hand, the oscillation period depends on 
the timescalcs of gene activation and deactivation. For 
example, the time t required to go from the fully acti- 
vated state Xi{t) ~ X* to Xi{t) ~ 1, where for a step-like 
repression function starts the target repression, is given 
by x*e^'^ = 1 — > r = ln{x*). This expression suggests 
that the promoter strength contributes logarithmically to 
the deactivation timescale, which is compatible with the 



dependence of the period on x* measured by numerical 
integration shown in Fig. [5J3. 



2. Increasing the growth rate can destabilize oscillations 

Given the promoter strength and the cooperativity 
of repression, defined by the specific properties of the 
genes and promoters implementing the repressilator, the 
growth-rate dependence of parameters defines a vertical 
path in the stability diagram of Fig. [2] (more precisely, the 
effective promoter strength F{pL) x* decreases nonlinearly 
with growth rate) . This path can cross the border of sta- 
bility of the limit cycle, defining a maximum growth rate 
at which stable oscillations can be substained. The pro- 
moter strength x* essentially defines the starting point 
of the path (more precisely the position at /i = 1 db/hr), 
while the function F(/i) is related to its length. Hence, 
both factors contribute to establish the range of growth 
rates in which oscillations are expected. 

Remarkably, the circuit can show qualitatively differ- 
ent dynamics depending on growth rate, and hence on 
nutrient conditions. Fig. U shows this for the case of 
cooperativity n = 3. Substained oscillations can be ob- 
served at slow growth, while in fast-growth conditions 
the dynamics can converge to a stable fixed point. The 
parameter range where oscillations are stable depends 
on where the repressilator is integrated through the gene 
dosage factor g appearing in F(/i) (sec Eq. [3]). 

If the repressilator is integrated on a plasmid, as the 
original in vivo experimental realization [l4|i the gene 
dosage g is simply given by the plasmid copy num- 
ber, which has a plasmid-speci fic g rowth-rate dependence 
that can be quite strong [l^, l29f . Fig. [J] compares the 
parameter regions corresponding to convergence to sta- 
ble oscillation or to a fixed point for a repressilator in- 
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FIG. 3: (Color online) Amplitude and period of oscilla- 
tion of a symmetric repressilator. The plots show os- 
cillation amplitude (A) and period (B) as a function of the 
promoter strength for three different values of cooperativity n. 
Both quantities are plotted in rescaled (dimensionless) units. 
The amplitude (concentration) is rescaled by the dissociation 
constant k, while the period has units of protein degradation 
rate /3p. For each value of n, Eqs. [7] were integrated numer- 
ically with values of a;* in the range 0.1 — 40 and step-size 
0.1. For each numerical solution, the oscillation amplitude 
and period were evaluated after convergence to a stable limit 
cycle. 



tegrated on the two plasmids Rl and pBR322 for which 
the copy number has been measured in different growth 
media ]W^ . 

The dosage of a chromosomal gene, instead, is deter- 
mined by its genomic position as set by the Cooper- 
Helmstetter model (see Appendix 1X1 and Fig.[Tj3). There- 
fore, for a repressilator integrated on the chromosome, 
the normalized circuit distance m from the replication 
origin (Ori) defines the growth-rate dependence of the 
dynamics (where m = represents Ori and m ~ I the 
replication terminus, Ter, Fig[Tj3). It should be noticed 
that the three genes composing the repressilator are as- 
sumed here to be inserted approximately in a single chro- 
mosomal location (or equivalently in different replichores, 
the oppositely replicated chromosome halves, but at the 
same distance m from Ori, see Fig.[T}3). The effects of the 
imbalance in gene dosage generated by different gene lo- 
cations on the chromosome are explored in scction llll B 21 
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FIG. 4: (Color online) The repressilator has different 
dynamics in different growth conditions. At fixed re- 
pression cooperativity (n = 3 in this plot), the repressilator 
can show stable oscillation or convergence to a steady state 
depending on the cell growth rate fi and the promoter strength 
X* (in units of k). The circuit substrate (plasmid, chromo- 
some) defines the range of growth rates at which a limit cycle 
can be observed. The cases of circuit integration on plas- 
mids {Rl and pBR322) and on chromosome near the origin 
of replication (m = 0) and between the origin and the termi- 
nus (m = 0.5) are represented in the plot. 



3. Period and amplitude of oscillation are growth rate 
dependent 

Fig. |4] shows that increasing the promoter strength a;* 
(which we remind also includes parameters not based on 
transcription) and the steepness of the repression func- 
tion n can make the limit cycle solution stable in a wide 
range of conditions. However, even if a repressilator is 
designed to exhibit oscillations in the experimentally ac- 
cessible conditions, the growth rate is still expected to 
influence the oscillation amplitude and period in a mea- 
surable way (i.e. by a factor that can exceed five, Fig.[5|). 
Indeed; the effective promoter strength F{fj,) x* decreases 
with increasing growth rate, because of the functional 
form of F{^), leading to a reduced amplitude and period 
of oscillation. Integrating Eqs. [7] numerically for parame- 
ter values corresponding to different growth rates allows 
to predict the oscillation period and amplitude for a chro- 
mosomally integrated repressilator or for a plasmid im- 
plementation, if the scaling of the plasmid copy number 
with growth rate is known. The presence of nonlincaritics 
in the system makes the dependence on growth rate of 
these two variables nontrivial, as represented in Fig. [5l 



B. The effect of intrinsic and position-induced 
asymmetry on the repressilator dynamics 

Differences in intrinsic properties of the genes compos- 
ing the circuit, such as affinity for RNA polymerase or 
ribosome binding, or gene dosage imbalances due to the 
specific gene location on the chromosome, lead to unequal 
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FIG. 5: (Color online) The oscillation amplitude and pe- 
riod of a symmetric repressilator vary significantly 
with growth rate. The plots show the scaling with growth 
rate of oscillation amplitude (concentration in units of k) and 
period (time in units of protein degradation rate) for a strong 
promoter (promoter strength x* = 70 in fc units) and coop- 
erativity n = 3. The three curves correspond to integration 
of the repressilator on plasmid Rl, plasmid pBR322, and to 
chromosomal integration of the three genes near the origin of 
replication (m = 0). A repressilator located near the repli- 
cation terminus (m = 1, not shown) follows a growth-rate 
dependence of amplitude and period similar to the case of 
plasmid pBR322. 



parameter values in the equations for the dynamics of Xi . 
This situation is generally referred to as "asymmetric". 
A certain degree of asymmetry is expected for generic 
synthetic realizations of the repressilator. However, with 
the exception of a few studies [1^, , its consequences 
on the dynamics have been poorly characterized theoret- 
ically Q. In the modeling framework adopted here, in- 
trinsic gene properties are recapitulated by the promoter 
strength. Additionally, using the Cooper-Helmstetter 
model (see Appendix]^ it is possible to account for the 
position-dependent scaling of gene dosage with growth 
rate. 

The two possible contributions to repressilator asym- 
metry, intrinsic gene properties and gene dosage, and 
their effects on the dynamics can be analyzed separately. 
We first address the dynamics at fixed growth rate of a 
repressilator composed of genes with different promoter 
strengths, and subsequently explore the consequences of 
position-induced asymmetry at different growth rates for 
a repressilator made of genes with equal intrinsic prop- 



erties, but different chromosomal location. 



1. Effects of asymmetry at fixed growth rate 

We consider the simplified situation in which only one 
of the genes of the repressilator differs from the others in 
its promoter strength by a factor A, but the three genes 
share the same growth rate dependence F{ijl). Thus, the 
single additional parameter A introduced in the model 
measures the level of circuit asymmetry. The system of 
equations describing the dynamics of a repressilator de- 
signed this way is 

Xi = X* F[y)R[x2) ~ xi 
xi = X* F{pL)R{xz) ~ xi 
X3 = A x*F{fi)R{xi)-X3 . (8) 

At fixed growth rate (for simplicity we take the case of 
growth rate fi — Idb/hr, where F(fj,) = 1), linear stability 
analysis can be applied to study the fixed point stabil- 
ity. As in the symmetric case described above, a Hopf 
bifurcation stands between the system convergence to a 
stable fixed point and the oscillatory solution. For each 
repression cooperativity level n, the stability diagram can 
be drawn as a function of the "basal" promoter strength 
X* and the asymmetry level A, as shown in Fig. [5] for 
n = 3. This diagram essentially shows that a high degree 
of asymmetry destabilizes the oscillations. Therefore, if 
the goal is to engineer a stable oscillator, a roughly sym- 
metric design is generally preferable. 

However, the minimum of the boundary curve between 
the two asymptotic dynamical behaviors (purple curve 
in Fig. [6]) docs not correspond to the symmetric case 
{Log{A) = 0). This result indicates that a symmetric 
system showing damped oscillations (as it is generally 
the case for parameter values just below the boundary in 
Fig. |6]) can be pushed toward a stable oscillation regime 
by slightly increasing the promoter strength of just one 
gene, if the resulting asymmetry is not too strong. 

Moreover, the presence of a gene with a different pro- 
moter strength breaks the symmetry in the dynamics of 
protein concentrations, making the oscillation amplitude 
(or the stable fixed point) gene-specific, as shown in Fig.[S] 
for two parameter sets. Thus, a simple test of the effec- 
tive symmetry of experimental repressilator realizations 
would entail measuring the oscillation amplitude of two 
fluorcsccntly-tagged protein products of genes in the cir- 
cuit. More specifically, the level of asymmetry, intro- 
duced by the presence of a gene with different intrinsic 
properties, affects the oscillation period and amplitude in 
a predictable way and with a gene-specific signature on 
the oscillation amplitude (Fig. [T]). Therefore, the global 
dynamical properties of the repressilator can be tuned 
simply by changing the parameters relative to a single 
gene. 
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FIG. 6: (Color online) Stability diagram of the asymmetric repressilator The parameter space is divided in the regions 
corresponding to convergence of the dynamics to a limit cycle and to a stable fixed point. The two parameters considered are 
the promoter strength x* (in units of k) and the level of asymmetry A (the factor by which the promoter strength of one of 
the genes differs from the two others). The cooperativity is fixed to n = 3 but the qualitative shape of the phase diagram 
holds for cooperativity values n > 2. While a high degree of asymmetry reduces the robustness of the oscillatory state, a slight 
asymmetry can stabilize the damped oscillations showed by the symmetric system. The right panel shows an example of the 
dynamics of the three protein concentrations for two sets of parameter values {{x* = 12, A — 10} and {x* = 4.5,^4 = 10}). 
The red continuous curve describes the dynamics of xi, the blue dashed curve refers to X2, and the orange dot-dashed to 
(same colors of the scheme in Fig. [TJ\). 
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FIG. 7: (Color online) Influence of asymmetry on pe- 
riod and amplitude of oscillation. The figure shows the 
nonlinear effects of the asymmetry induced by the presence 
of gene xs with a promoter strength that differs by a factor 
A from the other genes' promoter strength x* . The period 
(upper left plot) and the gene specific oscillation amplitudes 
(other plots) are shown as a function of the asymmetry level 
Log{A) for three promoter strength values that ensure oscil- 
lations in the all range of asymmetry explored. The curves 
are obtained by measuring period and amplitude of numerical 
solutions of Eq. [S]for values of Log{A) spaced by 0.01 in the 
range { — 1, 3}. 



The chromosomal position of genes affects the circuit 
dynamical properties at varying growth rates 



To isolate the effect of genomic position of genes on 
the dynamics of a chromosomally integrated repressila- 
tor, we consider the case of a circuit that is symmetric 
in terms of intrinsic gene properties, but where a gene is 
placed at a varying distance from the replication origin. 
This situation is illustrated in Fig. [5]4., where the genes 
with protein product concentrations xi and X2 are placed 
in contiguous positions, thus at approximately the same 
distance mi2 from Ori, while the gene corresponding to 
xg, is in position m^. In modeling terms, the equivalence 
of intrinsic gene properties corresponds to the assump- 
tion of identical effective promoter strengths XiFi{fi) at 
extremely slow growth, where the growth-dependence of 
gene dosage gi for each gene i is irrelevant, since gi — 1 
for every genomic position. In other words, the following 



relation is satisfied: x* F{fj, ~ 0) 



F3(/i ~ 0), where 



x" and F(fi) are the promoter strength and growth- 
dependence function of genes 1 and 2, and F3{fi) differs 
from F{p) only because of the position-dependent scaling 
of gene dosage with growth rate. 

It is straightforward to verify that the dynamics of a 
repressilator satisfying the above conditions is described 
by Eq. [8] with A = 2^'^'^™, where Am = mi2 - TO3 
and C is the time required for DNA replication (C ~ 
40 min for fast growth 3l|). The factor A is simply the 
ratio between the gene copy numbers, 53/51(2) as given 
by the Cooper- Helmstetter model (Appendix |^ , and 
encodes the growth-rate dependent level of asymmetry 
induced by gene position in an otherwise symmetrical 
circuit. Note that the sign of the relative position Am 
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FIG. 8: (Color online) Effect of gene-dosage induced asymmetry on oscillation stability (A) Scheme of the repressilator 
chromosomal configuration analyzed here. Two genes are placed at approximately the same distance mi2 from the origin of 
replication (possibly on opposite arms), while the third is in a different position 7713. This genomic configuration induces a 
growth-rate dependent copy-number imbalance between genes, as a consequence of DNA replication. The difference Am = 
mi2 — measures the effective asymmetry level (Am — corresponds to the symmetric case) since the ratio between the 
gene dosages (the factor A in Eq. [S]) is 2'^'^^'", where C is the DNA replication time. (B,C) The region of oscillation stability is 
plotted as a function of growth rate fi and relative position Am for different values of the promoter strength x*. The genomic 
configurations analyzed, and thus the values of Am, are obtained moving one gene along the chromosome while keeping the 
others near the replication origin (mi2 — 0) in plot B, or in a mid position (mi2 — 0.5) in plot C (see the corresponding schemes 
on the left). The plots show how gene positions infiuence the range of growth rates in which stable oscillations take place. (D,E) 
Effect of the relative position Am on the oscillation amplitude, normalized with its value in the symmetric case, of the three 
proteins is plotted for two different growth rates (promoter strength is x* = 75 in units of k). The sets of positions considered 
in D and E are the ones shown in plot B and C respectively. The change in oscillation amplitude that can be obtained by 
placing the same genes in different configurations depends on the growth conditions. 



has relevance. For example, the configuration with two 
genes near the replication origin and the third near the 
terminus (Am = — 1) leads to a quantitatively different 
dynamics with respect to the opposite configuration with 
Am = 1. The repressilator dynamics can be analyzed 
for different values of Am to explore the effects of gene 
position. 

Fig. [8j3 shows a stability diagram obtained moving 
along the chromosome one gene while the other two genes 
arc near the replication origin (mi2 = 0), for different lev- 
els of promoter strength x* . Analogously, the stability 
diagram in Fig. [8p is obtained keeping mi2 =0.5 and 
varying the third gene position m^. Interestingly, this 
analysis shows that at fixed growth rate the repressilator 
can converge to a limit cycle or to oscillations depend- 
ing on where genes are integreted on the chromosome, 
highlighting the importance of including the chromoso- 



mal gene coordinates as variables in genetic circuit mod- 
els. The positional effect is even more evident if two gene 
positions are varied (data not shown). 

As discussed in the previous section, the oscillation 
stability can be reinforced by increasing the promoter 
strength of one gene. A simple way to increase the ef- 
fective promoter strength of one gene for a chromosoma- 
Uy integrated circuit is moving its position towards Ori, 
since this increases its average gene copy number. Indeed, 
gene position has been used to modulate the oscillation 
features in an experimental synthetic implementation of 
an other genetic circuit [s^l- However, as we show for our 
case, the effect of the displacement of a gene is growth- 
rate dependent. In fact, while moving a gene from Ter 
to Ori allows an increment of its gene dosage of a fac- 
tor 4 for a /J, = 3 db/hr, this factor decreases with dou- 
bling time up to 1 for slow growth. The consequence of 
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this observation on the repressilator dynamics is evident 
looking at Fig. [5J3 and C. The same change in configu- 
ration (change in Am) has different effects depending on 
the growth rate. The modeling framework adopted here 
gives a quantitative prediction of the gain in oscillation 
stability that can be achieved by moving the genes along 
the chromosome depending on the experimental growth 
conditions, giving an experimental guideline for the best 
gene insertion sites in order to obtain the desired dynam- 
ics. 

Finally, Fig. [8jl) and E show the oscillation amplitude 
dependence on the relative gene positions for the config- 
urations described in Fig.|8j3 and C respectively. This de- 
pendence changes significantly at different growth rates, 
showing that the effect of gene position is strongly infiu- 
enccd by the cellular environment, and thus by growth 
rate. This feature is also relevant for synthetic biology, 
since identical experiments carried out with different nu- 
trient levels could in principle lead to different results. 



IV. DISCUSSION 

To sum up, this work addresses the dynamics of a 
paradigmatic bacterial oscillatory gene circuit, the re- 
pressilator, using for the first time a modeling frame- 
work that accounts for relevant physiological parameters, 
through growth-rate dependency [l^. From the mod- 
eling viewpoint, this framework entails assuming that 
the parameters of the dynamical system under exam are 
in fact dependent on a hidden "super-parameter" , the 
growth rate, which recapitulates cell physiology. Addi- 
tional models or experiments are required in order to 
obtain the dependence of each relevant parameter from 
the super-parameter. For the case of E. coU, all this in- 
formation is available. Since the parameters can move 
following the hidden variable, the phenomenology of the 
dynamical system with respect to the standard param- 
eters does not, by itself, allow to obtain the physically 
relevant phase diagram, if the behavior of the hidden 
variable is not known. This is important empirically, e.g. 
for comparing the results of experiment performed in dif- 
ferent conditions. 

A growth-rate induced "dynamical switch". Our results 
show that the dynamics is dependent on the growth rate 
in a fashion that is both qualitative and quantitative. 
Specifically, a symmetric repressilator will lose its oscilla- 
tory state with increasing growth rate unless its promoter 
strength is sufficiently high, and this phenomenology is 
expected to be observable in a wide range of experimen- 
tally accessible conditions. Indeed, the growth-mediated 
switch between different dynamics shown in Fig. |4] should 
be simple to observe experimentally, given the typical 
values of the parameters involved. The average protein 
copy numbers per cell span different orders of magnitude, 
from 10""'^ to 10^ [1^, while the values for the dissocia- 
tion constants have been reported to range between a 
few molecules [3J and a few thousands 341. With these 



numbers, the promoter strength values (i.e. protein con- 
centration at /i = 1 in units of k) analyzed here, such 
as X* G (0, 20) in the example in Fig. 21 are well in the 
physiological range. This suggests that a repressilator dy- 
namics characterized by loss of the oscillatory behavior at 
a critical growth rate should not be an exotic exception, 
and changes in both oscillation period and amplitude are 
likely to be measurable experimentally. Note however 
that some of the dissociation constant values reported 
in the literature were measured in vitro and thus could 
be underestimated with respect to in vivo values, where 
non-specific binding plays an important role. 

Substrate dependence of the dynamics. Further- 
more, the growth-rate dependent behavior of the cir- 
cuit is highly dependent on the substrate it is embed- 
ded in, a plasmid or a chromosome, and for a chromo- 
some on the detailed coordinates of the three promot- 
ers. As Fig. m and [S] suggest, in general the stability 
of oscillatory behavior is more sensitive to variations in 
growth rate for a repressilator integrated on plasmids 
than on the chromosome, as plasmid copy number can be 
strongly growth-rate dependent (as well as variable from 
cell to cell). However, integration on a high copy- number 
plasmid can naturally increase the protein concentration 
(thus the "promoter strength" parameter defined in this 
work), thus leading to more robust oscillations. There- 
fore, if the goal is engineering a stable synthetic oscillator, 
there is probably a trade-off between the advantage of an 
increased promoter strength typical of a high copy num- 
ber plasmid, and the unavoidable plasmid-specific growth 
rate dependence (and cell-to-cell variability) of the copy 
number. Accounting for these variables in a synthetic 
biology context could be rather difficult, since measure- 
ments of plasmid average copy number and in particular 
cell-to-cell variations are not easy [s^ , and there exists a 
wide range of behaviors in plasmids, depending on their 
biological features and the mechanism that they use for 
segregation. 

On the other hand, for chromosomally integrated re- 
pressilators, the dynamics depends nontrivially on gene 
position. Recently, it has been shown that the spatial 
ordering of a set of "important" genes along the chro- 
mosome is strongly conserved between different bacte- 
rial species and largely corresponds to their expression 
pattern during growth [36j . pointing to a functional role 
for gene chromosomal position. The example of a puta- 
tive chromosome-integrated repressilator analyzed here 
suggests that the dynamics of genetic networks, in fast- 
growing bacteria, should be influenced by the genomic 
position of its components. In fact, circuit behavior is 
affected by growth rate in a position-specific manner, be- 
cause of gene dosage effects created by overlapping repli- 
cation rounds. More generally, it's tempting to spec- 
ulate that the evolutionary pressure to keep a certain 
gene order with respect to genome replication may be 
partially due to natural selection of specific network dy- 
namics defined by the combination of gene positions and 
cell growth state. 
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In this perspective, the analysis of the phenotypic con- 
sequences of chromosomal rearrangements, such as large 
inversions [s^], should be revisited taking into account 
the growth conditions. In fact, as shown in Fig.[5]for the 
repressilator case, a variation in the relative position of 
genes involved in a regulatory circuit can have different 
consequences depending on the growth rate. 

Role of noise. Some relevant considerations can be 
made about the possible role of noise in the results given 
here within a purely deterministic framework. In general, 
noise can strongly affect the dynamics of a repressilator. 
For example, in the in vivo realization of the repressila- 
tor [iJl only about 40% of the cells displayed oscillations, 
with high cell-to-cell variability in oscillation period and 
amplitude. Several studies have analyzed the possible 
impact of noise, focusing on the stochasticity that arises 
from the discrete nature of the molecular players and 
from the inherent randomness of their interactions (to- 
gether refereed to as intrinsic noise) [l3, [IB, 113 • The 
main result is that intrinsic noise can both play a con- 
structive role in oscillation robustness and a destruc- 
tive one. The constructive phenomenon can enlarge the 
parameter space of oscillations through a resonance ef- 
fect [13 . The destructive one causes strong cell-to-cell 
variability in oscillation amplitude and period [l^ . 

While intrinsic noise can be a relevant factor and could 
partially explain the experimentally observed variability, 
the dominant source of noise might be due to fluctua- 
tions in global cellular parameters (extrinsic noise) ^ such 
as ribosome or polymerases concentration. This has been 
shown to be the case in E. coli, for relatively high expres- 
sion (more than approximately 10 proteins per cell) [23j . 
Since oscillations in the repressilator generally require 
strong promoters, the variability in the circuit dynamics 
is expected to be highly sensitive to the extrinsic noise 
level, which adds up, for plasmids, to the aforementioned 
cell-to-cell copy-number variation. These considerations 
point to an interest in considering the stochastic aspects 
of the circuit. However, in order to extend the mean- 
field model introduced here to analyze the growth rate 
dependence of the cell-to-cell variability of the repressi- 
lator dynamics, it would be necessary to know how the 
extrinsic noise scales with growth rate. Unfortunately, 
there are no experimental data concerning this scaling, 
making the extension of this work to the stochastic case 
premature. 

Nevertheless, the possibility of a switch between dif- 
ferent dynamical regimes in response to the physiologi- 
cal cell state opens interesting considerations about the 
robustness of this oscillatory genetic circuit. Fluctua- 
tions in physiological parameters such as the growth rate 
fall in the broadly-defined category of extrinsic noise. It 
has been suggested that fluctuations in the growth rate, 
mainly through its influence on protein dilution, can ac- 
count for a considerable part of the measured extrinsic 
noise [H, . Our analysis suggests that fluctuations in 
cell parameters linked to growth can introduce cell-to-cell 
variability in the circuit dynamics (convergence to oscil- 



lation or to a stable fixed point) as well as in oscillation 
period and amplitude. Therefore, the reasons behind the 
lack of robustness of the repressilator realized in vivo [U 
should be also searched in the varibility of physiological 
parameters rather than focusing exclusively on intrinsic 
noise effects. 

Connection to experiments and outlook. Finally, we 
believe these findings could be relevant from both a sys- 
tems biology and a synthetic biology perspective. There 
is a long list of endogenous oscillators in bacteria [i^ 
that are interesting for the former discipline, and need to 
be understood within the framework adopted here. The 
most important examples are circadian clocks and the 
cell cycle itself. We previously studied the dynamics of 
the DnaA oscillatory circuit, which is determinant in this 
last process [l^ . In this case, the timescale of oscillations 
matches (by definition) the cell cycle time, thus the ap- 
proximations defined in section fll Al are not valid. More 
complex models arc required, and there is no commonly 
accepted reference for this description. However, it is 
interesting to note that in the simple modeling frame- 
work adopted here, the period of oscillation decreases 
with growth rate without the need of specific additional 
regulation. In absence of overlapping replication rounds, 
this is exactly the kind of behavior behavior desired for 
an oscillator regulating the triggering of DNA replication, 
such as the DnaA circuit: a shortening of the initiation 
time is required when the cell volume grows faster to 
synchronize DNA replication and volume doubling (the 
situation becomes more complex at fast growth [22l|'). 

In contrast, circadian clocks need to be resilient to 
changes in the cell doubling time, and thus in the growth 
rate, in order to keep a steady 24-hour period in vari- 
able environmental conditions, and thus cannot measure 
time using the cell cycle. The consequent decoupling be- 
tween the cell cycle and the circadian rhythm has in- 
deed been verified in cyanobacteria [4l| - |43j . Our results 
suggest that the dynamics of a genetic oscillator is nat- 
urally strongly connected to the cellular growth rate. 
Therefore, specific regulatory mechanisms are required 
in order to compensate for these effects and render a 
circadian oscillator insulated from the growth state. Al- 
though circadian clocks appear to be primarily based on 
post-translational circuits in bacteria (40| . the proteins 
involved are the result of a gene expression process, and 
thus in principle coupled with growth rate [15| . It would 
be interesting to evaluate experimentally if the promoters 
regulating these proteins are more buffered as a function 
of growth rate compared to others. Quantitative models 
taking into account the parameter dependence on growth 
rate, such as the one presented here, may be important 
to pose the question of the growth-rate robustness of the 
circadian cycle. For example, the circuit architectures 
and the type of regulations selected by evolution to com- 
pose circadian oscillators might be, at least in part, con- 
strained by the implementation of the observed growth- 
rate independence. 

Finally, from a synthetic biology standpoint, changing 
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the conditions in which the cells are grown alters quanti- 
tatively the characteristics of the repressilator dynamics 
in a predictable manner. This offers the possibility of 
external control of the circuit behavior by simply oper- 
ating on macroscopic variable related to physiology such 
as the type of nutrient supply or the temperature. This 
way, the engineering and control of the dynamics can 
be performed by tuning environmental conditions in a 
model-guided way, rather then by modifing the genetic 
components, which can be technically complex. More- 
over, the modeling strategy presented here can suggest 
the best chromosomal position for gene insertion depend- 
ing on the specificities of the desired circuit dynamics. 




^ ^ 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 

growth rate |J [db/hr] 
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FIG. 9: (Color online) Cell cycle and population average 
of gene dosage. The gene dosage g averaged over the cell 
cycle (dashed blue lines) and averaged over a population in 
balanced exponential growth (continuous red lines) are shown 
as a function of the growth rate for different chromosomal po- 
sitions m (normalized distance from Ori). In the physiological 
range of growth rates, the two quantities do not present a sub- 
stantial quantitative difference. 



1. Cooper-Helmstetter model and gene dosage 



Appendix A: Population averages vs cell-cycle 
averages 

A common assumption in genetic circuit modeling is 
that the contribution of the cell-cycle dynamics can be 
neglected, at least when one is interested in time scales 
longer than the doubling time, which is often the case 
given the typical high protein stability. This approxi- 
mation allows the use of effective parameters obtained 
averaging over the cell cycle [l^l- However, as discussed 
in the text, the growth-rate dependences of gene expres- 
sion parameters are derived from experimental observa- 
tions of average cellular properties in a population [ig} . 
which are also affected by the cell age distribution. This 
Appendix shows the quantitative difference between pop- 
ulation and cell-cycle averages for two quantities that are 
well characterized, the gene dosage due to DNA replica- 
tion and the cellular volume. We will show that in these 
two cases population averages can be used in dynamic 
models of genetic circuits without introducing significant 
errors. The assumption that this result can be gener- 
alized to other quantities justifies, although not rigor- 
ously, the use of available experimentally estimated pop- 
ulation averages in dynamic models for those quantities 
whose time dependence (or even cell-cycle averages) are 
not known [l^. This is the case for ribosome or poly- 
merase concentrations at different growth rates, which 
are crucial to determine the growth-rate dependence of 
transcription and translation rates in Eq. |31 



DNA replication in fast-growing bacteria such as 
E. coli typically starts from a single replication ori- 
gin (Ori) and proceeds bidirectionally along the circu- 
lar chromosome until it reaches the replication terminus 
(Ter). The Cooper-Helmstetter model [2l| establishes 
the relation between growth-rate and replication timing 
such that DNA copies are produced on time for each new- 
born cell. The model is based on the empirical observa- 
tion that the time necessary for chromosome replication 
(called "C period") and the time period between com- 
pletion of chromosome replication and the following cell 
division (D period) are approximately constant (at least 
for fast-dividing cells [Slj, with doubling times less than 
Ihr). Since at time C + D the cell divides, a time lag X 
before initiation is necessary to make the total replica- 
tion time X + C + D an integer multiple of the doubling 
time r, "synchronizing" DNA replication and cell divi- 
sion. Thus, the following relation has to be satisfied 

X + C + D = {n+l)T , (Al) 

where n = Int[ '-'~^^ ] is the integer number of times that t 
divides C+D. Starting from this relation, it can be easily 
shown that the number of origins present at initiation is 
exactly 2" [l^, [2l|. More generally, we can consider a 
gene at a chromosomal position defined by its normalized 
distance m from Ori, i.e. m = represents a gene in Ori 
and m 1 in Ter. The copy number of this gene, g, 
changes during the cell cycle following 

. s ^ / 2"' if < t < (n' + l)r - (C(l - m) + D) 
' \ 2"'+i if {n' + 1)t - (C(l - rn) + D) < t < t 
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where 



Int 



C(l-m)+_Di 



Therefore, the gene dosage 



averaged over the cell cycle (which could be measured 
following a single cell lineage and averaging over time), 
is given by 



(git)) cell cycle 
= 2"' [1 



g{t)dt 



1 

T JO 

n' + fi{C{l-m) + D)] 



(A2) 



On the other hand, the population age structure must 
be taken into account when evaluating the average gene 
dosage in a cell population. For ideal "balanced expo- 
nential" growth with rate fi this distribution is given 

by Mm 



Thus, the population average is 



(A3) 



(git)) 



population 



ait, fi)g{t)dt ~ 



^ 2t^lC(l-m)+D] 



(A4) 



This is the expression typically used to evaluate the gene 
dosage [H, [la|. As shown in Fig. ^ the difference be- 
tween cell-cycle averages (dashed blue lines) and popula- 



tion averages (red continuous lines) in the physiological 
range of growth rates is negligible for all gene positions. 



2. Cell volume growth 

Similar considerations can be carried out in the case 
of the average cell volume. The functional form of the 
volume increase in time during a cell cycle has long 
been debated [1^ , with two prevailing hypotheses of lin- 
ear growth (constant rate) or exponential growth (size- 
dependent rate), although more complex dependences 
have been proposed [1^. Recent experiments strongly 
suggest an exponential growth |47| , and we assumed this 
functional form (note that the same reasoning could be 
applied to linear growth straightforwardly, so this choice 
has no consequences on any of the results). With a vol- 
ume growth of the form V{t) = Vb2^*, the mean volume 
over a cell cycle is {V(t))ceU cycle = Vb/Zn2, while the in- 
tegration over the population leads to {V {t)) population = 
2 ln2 Vq. All the volume growth-rate dependence is hid- 
den in Vq, and experimental results indicate that this 
dependence is approximately exponential j48j . The dif- 
ferent numerical factors introduced by averaging over the 
cell cycle or over the population do not affect significantly 
the volume growth-rate dependence. 



[1] A. Goldbeter, Biochemical Oscillations and Cellular 
Rhythms (Cambridge University Press, 1997). 

[2] A. Goldbeter, Nature 420, 238 (2 002), URL 
|http : //dx ■ doi . org/10 . 1 038/nature01259 

[3] J. J. Tyson, R. Albert, A. Goldbeter, R Ruoff, and 
J. Sible, J R Soc Interface 5 Suppl 1, SI (2008), URL 
http : //dx . doi . org/10 . 1098/rsif . 2008 .0179 .focus 

[4] O. Purcell, NT^jTSavery, CTS. Grierson, and M. di 
Bernardo, J R Soc Interface 7, 1503 (201 0), URL 
|http : //dx ■ doi ■ org/10 ■ 1098/rsif .2010.01831 

[5] L. Bintu, N. E. Buchler, H. G. Garcia, U. Ger- 
land, T. Hwa, J. Kondev, and R. Phillips, 
Curr Opin Genet Dev 15, 116 (2005), URL 
,http://dx . doi . org/10 . 1016/ j . gde . 2005 . 02 . 007 

[6] U. Alon, An Introduction to Systems Biology: Design 
Principles of Biological Circuits (Chapman & Hall/CRC 
Mathematical & Computational Biology, 2006). 

[7] M. E. Wall, W. S. Hlavacek, and M. A. Sav- 
ageau, Nat Rev Genet 5, 34 (2004), URL 
|http : //dx ■ doi . org/10 ■ 1038/nrgl244| 

[8] G.-W. Li and X. S. Xie, Nature 475, 308 (2011), URL 
|http : //dx ■ doi . org/10 . 1038/naturel0315 i 

[9] E. Andrianantoandro, S. Basu, D. K. Karig, and 
R. Weiss, Mol Syst Biol 2, 2006.0028 (2006), URL 
http : //dx . doi . org/10 . 1038/msb4100073 
[10] N. J. Guido, X. Wang, D. Adalstcinsson, D. McMillen, 
J. Hasty, C. R. Cantor, T. C. Elston, and 
J. J. Collins, Nature 439, 856 (2 006), URL 
|http : //dx ■ d o i . org/10 . 1038/nature04473 | 
[11] J. Haseloff and J. Ajioka, J R Soc In- 
terface 6 Suppl 4, S389 (2009), URL 



|http : //dx ■ doi ■ org/10 ■ 1098/rsif . 2009 .0176 ■ focus [ 
[12] N. Nandagopal and M. B. Elowitz, 

Science 333, 1244 (2011), URL 

|http : //dx ■ doi ■ org/10 ■ 1 126/science . 1207084] 

[13] W. Weber and M. Fussenegger, Nat Rev Genet 13, 21 
(2012) , URL http : //dx . doi . org/10 . 1038/nrg3094 

[14] M. B . Elowitz and S. Leibler, Nature 403, 335 (2000), 
URL [http : //dx ■ doi ■ org/10 .1038/35002125 . 

[15] S. Klumpp, Z. Zhang, and T. Hwa, Cell 139, 1366 (2009). 

[16] H. Bremer and P. Dennis, In: F.C. Neidhardt, Editor, 
Escherichia coli and Salmonella, ASM Press, Washing- 
ton, D.C. pp. 1553-1569. (1996). 

[17] M. Scott, C. W. Gunderson, E. M. Mateescu, Z. Zhang, 
and T. Hwa, Science 330, 1099 (2010). 

[18] M. Scott and T. Hwa, Curr Opin Biotechnol (2011). 

[19] C. Tan, P. Marguet, and L. You, 
Nal _ Chem Biol 5, 84 2 _ (2009 ), URL 
htt^pT7/dx. doi.org/107l638/nchemb 

[20] R. Kwok, Nature 463, 288 (2010), URL 
[http : //dx . doi . org/lb ■ i0 38/463288a 

[21] S. Cooper and C. Helmstetter, Journal of Molecular Bi- 
ology 31, 519 (1968). 

[22] M. A. A. Grant, C. Saggioro, U. Ferrari, 
B. Bassetti, B. Sclavi, and M. C. Lago- 
marsino, BMC Syst Biol 5, 201 (2011), URL 
http : //dx ■ doi . org/10 . 1186/1752-0509-5-20 1 

[23] Taniguchi, R I Cim, G7-W. Li, 

H. Chen, M. Babu, J. Hearn, A. Emili, and 
X. S. Xie, Science 329, 533 (2010), URL 
http : //dx . doi . org/10 . 1126/science . 1188368 

[24] R. Marathe, V. Bierbaum, D. Gomez, and S. Klumpp, 



13 



J Stat Phys pp. 1-20 (2012), ISSN 0022-4715, [37 
10.1007/sl0955-012-0459-0. 

S. MiiUer, J. Hofbauer, L. Endler, C. Flamm, S. Wid- 
der, and P. Schuster, J Math Biol 53, 905 (2006), URL [38 
|http : //dx ■ doi . org/10 . 1007/s00285-006-0035-9| 
N. Strelkowa and M. Barahona, J R 
Soc Interface 7, 1071 (2010), URL [39 
|http : //dx ■ doi . org/lo . 10987rsif . 2009 . 0487 1 
O. Buse, R. Prez, and A. Kuznetsov, Phys Rev E Stat 
Nonlin Soft Matter Phys 81, 066206 (2010). 
R. Blossey and C. V. Giuraniuc, Phys Rev E Stat Nonhn [40 
Soft Matter Phys 78, 031909 (2008). 
S. Klumpp, PLoS One 6, e20403 (2011), URL 
http : //dx ■ do i ■ org/10 . 1371/ j ournal . pone . 0020403 [41 
H. El Samad, D. Del Vecchio, and M. Khammash, in 
Proceedings of the 2005 American Control Conference [42 
(ieee Cat. No. 05ch37668) (American Automatic Control 
Council, 2005). 

O. Michelsen, M. J. T. de Mattos, P. R. Jensen, and F. G. [43 
Hansen, Microbiology 149, 1001 (2003). 
M. R. Atkinson, M. A. Savageau, J. T. Myers, and A. J. 
Ninfa, Cell 113, 597 (2003). [44 
L. Bintu, N. E. Buchler, H. G. Garcia, U. Ger- [45^ 
land, T. Hwa, J. Kondev, T. Kuhlman, and 
R. Phillips, Curr Opin Genet Dev 15, 125 (2005), URL [46 
http://dx. doi.org/10. 1016/ j .gde . 2005 . 02 . 006 
Y. Setty, A. E. Mayo, M. G. Surette, and U7 Alon, 
Proc Natl Acad Sci U S A 100, 7702 (2003 ), URL [47 
|http : //dx ■ doi ■ org/10 . 1073/pnas . 1230759 100 | 
S. Tal and J. Paulsson, Plasmid' 67, l67 (2012 ), URL 
http://dx.doi.org/10. 1016/ j ■plasmid.2012. 01 .004| 
P. Sobetzko, A. Travers, and G. Muskhelishvili, Proc [ 
Natl Acad Sci U S A 109, E42 (2012) URL 
|http : //dx . doi . org/10 . 1073/pnas . 1108229109| 



E. Esnault, M. Valens, O. EspeU, and F. Boc- 
card, PLoS Genet 3, e226 (2007), UR L 
|http : //dx . doi ■ org/ 10 ■ 137 1/ j ournal . pgen . 0030226 [ 
M. Yoda, T. Ushikubo, W. Inoue, and M. Sa- 
sai, J Chem Phys 126, 115J101 (2007), URL 
'http : //dx . doi . org/lO . 1063/172539537 
S. Tsuru, J. Ichinose, A. Kashiwagi, 



W. Ying, 
Phys Biol 



K. 



Kaneko, 
036015 



B.- 

and T. Yomo, 
(2009), URL 



|http://dx.doi. org/10 .1088/1478-3975/6/3/0360151 

P. Lenz and L. S0gaard-Andersen, Nat 

Rev Microbiol 9, 565 (2011), URL 

http : //dx . doi . org/10 . 1038/nrmicro2612 

T. Kondo, T. Mori, N. V. Lebedeva, S. Aoki, M. Ishiura, 

and S. S. Golden, Science 275, 224 (1997). 

T. Mori and C. H. Johnson, J 

Bacteriol 183, 2439 (2001), URL 

|http : / /dx ■ doi ■ org/10 . 1 128/ JB . 183 . 8 . 2439- 2444 . 2001 [ 



I. Mihalcescu, W. Hsing, 

S. Leibler, Nature 430, 81 (2004), 
|http : //dx . doi . org/10 ■ 103 8/nature02533 
E. O. Powell, J Gen Microbiol 15, 492 (1956). 
A. L. Koch, Grit Rev M icrobiol 1 9, 17 (19 93), 
|http://dx. doi. org/ 10731 09/10408 419309113 521 
G. Reshes, S. Vanounou, I. Fishov, 
M. Feingold, Biophys J 94, 251 (2008), 



and 
URL 



URL 

and 
URL 



|http : //dx . doi . org/10 ■ 1529/biophys j . 107 . 104398 | 
M. Godin, F. F. Delgado, S. Son, W. H. Grover, A. K. 
Bryan, A. Tzur, P. Jorgensen, K. Payer, A. D. Grossman, 
M. W. Kirschner, et al., Nat Methods 7, 387 (2010), URL 
http : //dx . doi . org/10 . 1038/nmeth . 1452 
M. Schaechter, O. Maaloe, and N. O. Kjeldgaard, J Gen 
Microbiol 19, 592 (1958). 



